GC-MS/MS Method for Determination of Polycyclic Aromatic Hydrocarbons in Herbal Medicines

Polycyclic aromatic hydrocarbons (PAHs) are hydrophobic organic contaminants that have a highly carcinogenic and mutagenic nature. This study aimed to develop and validate a sensitive analytical method to determine 8 PAHs in 51 herbal medicines (HMs) using gas chromatography (GC)-tandem mass spectrometry (MS/MS). Liquid––liquid extraction and florisil SPE cartridge purification were basically adopted for pretreatment. For the samples containing essential oil, starch grain, etc., N,N-dimethyl formamide/water mixture (9:1, v/v) was added in the extraction step. The multiple reaction monitoring (MRM) conditions were newly obtained by the infusion of reference solutions of the targeted compounds at a concentration of 100 ng/mL into the GC-MS/MS system used in this study. The 51 items were classified according to whether or not they contained essential oil. Eight PAHs were not detected in 39 (8.3%) of the 459 samples monitored. The total content of 8 PAHs ranged from 0.45 μg/kg in Anemarrhenae Rhizoma to 270.94 μg/kg in Zingiberis Rhizoma. The average content of those ranged from 0.9 μg/kg in Araliae Continentalis Radix to 110.8 μg/kg in Coptidis Rhizoma Preparata cum Vinum. The results of this study prove that the proposed method is useful for determining 8 PAHs in HMs.


Introduction
Polycyclic aromatic hydrocarbons (PAHs) are hydrophobic organic contaminants that are composed of two or more fused aromatic rings [1]. Once formed by incomplete combustion, PAHs spread to all environmental components such as air, water, sediment, soil and plants, polluting aquatic and terrestrial species [2]. In 1984, the United States Environmental Protection Agency (USEPA) designated 16 PAHs as priority environmental pollutants of carcinogenic or mutagenic [3]. In 2007, the European Food Safety Authority (EFSA) reported descriptive statistics for the concentrations (µg/kg) of up to 16 PAHs in 9714 food products [4]. According to the report, 3086 samples (31.8%) did not contain any PAH above the limit of detection. In 2008, the EFSA concluded that benzo[a]pyrene (BaP) is not a suitable indicator for the presence of PAHs in food and that both 8 PAHs (benzo[a]anthracene (BaA), chrysene(CHR), benzo[b]fluoranthene (BbF), benzo[k]fluoranthene (BkF), BaP, dibenzo[a,h]anthracene (DahA), benzo[g,h,i]perylene (BghiP) and indeno [1,2,3-c,d]pyrene (IcdP)) and 4 PAHs (BaA, CHR, BbF and BaP) are the most suitable markers of PAHs in food [5][6][7]. The International Agency for Research on Cancer (IARC) classified BaP as Group 1, DahA as Group 2A, BaA, CHR, BbF, BkF and IcdP as Group 2B and BghiP as Group 3 [8].
Since ancient times, herbal medicines (HMs) have been used for various purposes such as food ingredients or therapeutic treatments. In the modern era, although synthetic drugs are effective for various diseases, the interest in substances with the therapeutic effects of plants grown in nature is increasing. However, HMs are facing stumbling in the case of some samples containing essential oils, there were some items that had good peaks and recoveries without interfering substances even when extracted only with hexane (non-DFA extraction method), and the DFA extraction method was not required for them. Meanwhile, the MRM conditions for 8 PAHs were set as compiled in Table 1 and were slightly different from those in the study of Anderson et al. [34] in terms of the m/z values. The determined MRM mode was found to have high sensitivity for the 8 PAHs in the GC-MS/MS analysis. Consequently, this study set up well-suited quantitative ions and confirmation ions in the MRM mode for the identification and quantification of PAHs in the HMs.  Sanguisorbae Radix Carbonisatum, which showed stable chromatogram baseline and relative low detected concentrations of 8 PAHs, was selected as a representative of the first group of 19 thermal processed HMs. Moreover, all of the first groups showed good peaks in the circumstance of the non-DFA extraction method. The second group, consisting of 13 natural HMs, was classified according to whether or not they contained essential oil, starch grain and calcium oxalate by referring to the properties and characteristics described in the Korean pharmacopoeia. Although Puerariae Radix, Glycyrrhizae Radix et Rhizoma, Atractylodis Rhizoma Alba, Bupleuri Radix, Achyranthis Radix, Ginseng Radix, Cnidii Rhizoma and Cyperi Rhizoma contained essential oil and starch grain, they showed good peaks even by the non-DFA extraction method. So, only Angelicae Gigantis Radix, which has starch granules in its oil cell, was chosen as a separate representative to which the DFA extraction method was applied. An important reason for this choice was that, for instance, Angelicae Gigantis Radix could not be ionized when pretreated by using the non-DFA extraction method ( Figure 2). By comparing the blanks of the second group, Achyranthis Radix and Scutellariae Radix, along with Angelicae Gigantis Radix, were all selected as representatives because they exhibited stable chromatogram baselines and relatively low detected concentrations of 8 PAHs. For the third group, consisting of 19 natural HMs, the DFA extraction method was applied for Asiasari Radix et Rhizoma, Acori Graminei Rhizoma, Atractylodis Rhizoma and Araliae Continentalis Radix, which contain essential oil because of the properties and characteristics described in the Korean pharmacopoeia. For example, Araliae Continentalis Radix showed a poor peak for CHR-12 when extracted with the non-DFA extraction method but exhibited a good peak after applying the DFA extraction method ( Figure 3). By comparing the blanks of the third group, Asiasari Radix et Rhizoma and Anemarrhenae Rhizoma among 19 HMs of the third group showing stable baselines and relatively low detected concentrations of 8 PAHs, were selected as representatives. Consequently, as for six representatives, the validation factors such as selectivity, linearity, LOD, LOQ and recovery rates were examined and illustrated in following sections.

Matrix Effect
Depending on the matrix of the blank sample, the ionization efficiency can be lower or higher while mass spectrometry is being used [35]. Since most HMs originate from natural plants containing various compounds, which may affect the quantitative measurements, matrix effects are inevitable when analyzing substances or contaminants they contain. Therefore, there is a strong matrix effect that causes a difference in ionization efficiency. To reduce the matrix effect, a calibration curve is prepared using the matrixmatched method in which a standard product is dissolved in a blank sample. That is, when there is a matrix effect, linearity should be evaluated using a standard solution containing the sample matrix [36]. Matrix effect can be classified as strong matrix effect below 50% or above 150%, moderate matrix effect between 50-80% or 120-150%, and low matrix effect between 80-120% [37]. Taking the criteria into account, the matrix effect on 8 PAHs was investigated. The matrix effects for the six representatives are summarized in Table  2. In Sanguisorbae Radix Carbonisatum, although IcdP and DahA had strong matrix effects of 189% and 167%, respectively, those showed acceptable recoveries. Except for IcdP and DahA of Aconiti Ciliare Tuber Preparata and Morindae Radix Preparata cum Vinum, most of the other items showed low to moderate matrix effects on all PAHs (Table S1). It is also noteworthy that, owing to acceptable recoveries that were shown for most of the items, calibration curves could be determined to be prepared by dissolving a standard internal material in a standard solution instead of using the matrix-matched method.

Matrix Effect
Depending on the matrix of the blank sample, the ionization efficiency can be lower or higher while mass spectrometry is being used [35]. Since most HMs originate from natural plants containing various compounds, which may affect the quantitative measurements, matrix effects are inevitable when analyzing substances or contaminants they contain. Therefore, there is a strong matrix effect that causes a difference in ionization efficiency. To reduce the matrix effect, a calibration curve is prepared using the matrix-matched method in which a standard product is dissolved in a blank sample. That is, when there is a matrix effect, linearity should be evaluated using a standard solution containing the sample matrix [36]. Matrix effect can be classified as strong matrix effect below 50% or above 150%, moderate matrix effect between 50-80% or 120-150%, and low matrix effect between 80-120% [37]. Taking the criteria into account, the matrix effect on 8 PAHs was investigated. The matrix effects for the six representatives are summarized in Table 2. In Sanguisorbae Radix Carbonisatum, although IcdP and DahA had strong matrix effects of 189% and 167%, respectively, those showed acceptable recoveries. Except for IcdP and DahA of Aconiti Ciliare Tuber Preparata and Morindae Radix Preparata cum Vinum, most of the other items showed low to moderate matrix effects on all PAHs (Table S1). It is also noteworthy that, owing to acceptable recoveries that were shown for most of the items, calibration curves could be determined to be prepared by dissolving a standard internal material in a standard solution instead of using the matrix-matched method.

Linearity and Sensitivity
The linearity between the analytical concentration and the chromatogram peak was evaluated by the correlation coefficient (R 2 ). The R 2 values for 8 PAHs in the quantitative range (1-100 ng/mL) were greater than 0.998 (Table 1), indicating that the analytical method proposed in this study provided good linearity due to its good selectivity (Figure 4). LOD and LOQ for the 8 PAHs in each of the six representatives were summarized in Table 2 122.2%, with the intra-day RSD of 1.5-19.0% and the inter-day RSDr of 1.0-25.1%. Inter-205 estingly, both the minimum and maximum values of recovery rates were found in Scutel-206 lariae Radix among the six representatives. In detail, the highest recovery rate was 125.0% 207 at 2 LOQ for BaA, and the lowest recovery rate was 54.7% at 10 LOQ for CHR (Table S2). 208 It is speculated that this is probably due to the characteristic matrix of Scutellariae Radix. 209 In the near future, analytical studies on the differences in matrices of HMs by nutritional 210 analysts are required.  As stated above, the validation factors obtained through these verified methods con-219 formed to the AOAC guidelines [38]. Particularly, the recoveries of PAHs at low, medium 220 and high concentrations appeared to be within an acceptable range, although Scutellariae 221 Radix, Anemarrhenae Rhizoma and Asiasari Radix et Rhizoma showed slightly overesti-222 mated recovery values. In addition, the RSDs were entirely reasonable and satisfactory 223

Recovery and Precision
Standard substances corresponding to three different concentrations (low, medium, and high), such as 2, 5, and 10 times the LOQ, were added to the samples within the quantification range of six representatives. The accuracy of 8 PAHs was evaluated with the recovery rate calculated from the measured results. The precision was evaluated as the relative standard deviation (RSD) of the measured results. The recovery rates and RSDs are described in Table 2. As for Sanguisorbae Radix Carbonisatum, the recovery rate of PAHs was 66.4-104.7%, with the intra-day RSD of 0.8-22.3% and the inter-day RSDr of 1.0-17.1%. The recovery rate from Achyranthis Radix was 71.2-106.5%, with the intra-day RSD of 0.5-9.8% and the inter-day RSDr of 0.5-15.1%. The recovery rate from Scutellariae Radix was 54.7-125.0%, with the intra-day RSD of 0.4-12.3% and the inter-day RSDr of 3.1-32.1%. The recovery rate from Angelicae Gigantis Radix was 70.7-94.3%, with the intra-day RSD of 0.6-14.3% and the inter-day RSDr of 1.1-22.9%. The recovery rate from Anemarrhenae Rhizoma was 83.7-124.4%, with the intra-day RSD of 1.2-13.7% and the inter-day RSDr of 0.3-17.2%. The recovery rate from Asiasari Radix et Rhizoma was 70.2-122.2%, with the intra-day RSD of 1.5-19.0% and the inter-day RSDr of 1.0-25.1%. Interestingly, both the minimum and maximum values of recovery rates were found in Scutellariae Radix among the six representatives. In detail, the highest recovery rate was 125.0% at 2 LOQ for BaA, and the lowest recovery rate was 54.7% at 10 LOQ for CHR (Table S2). It is speculated that this is probably due to the characteristic matrix of Scutellariae Radix. In the near future, analytical studies on the differences in matrices of HMs by nutritional analysts are required.
As stated above, the validation factors obtained through these verified methods conformed to the AOAC guidelines [38]. Particularly, the recoveries of PAHs at low, medium and high concentrations appeared to be within an acceptable range, although Scutellariae Radix, Anemarrhenae Rhizoma and Asiasari Radix et Rhizoma showed slightly overestimated recovery values. In addition, the RSDs were entirely reasonable and satisfactory values.

Cross Validation
A cross-validation experiment was performed to validate the proposed analytical method by comparatively evaluating inter-laboratory selectivity, linearity, LOD, LOQ, recovery and precision. For this purpose, three different concentrations (low, medium and high concentrations set as 2, 5 and 10 times the LOQ, respectively) of standard substances were added to the six representatives. As a result, no disturbed peaks appeared in the chromatogram of GC-MS/MS. In all cases, the linearity was greater than 0.994, and the LOQs for the 8 PAHs in each of the six representatives were in the range of 0.44 and 1.18 µg/kg ( Table 3) (Table 3). Since the recoveries and RSDs of most samples obtained through cross-validation satisfied the acceptance criteria of AOAC [38], the proposed analytical method was considered to be useful for determining 8 PAHs in six representatives.

Applicability
This newly developed method was further tested by applying it to a total of 45 items, excluding the six representatives analyzed during the validation process presented in this study. The recoveries were tested at two times the LOQ, whereas Linderae Radix (BbF, BkF, BaP) and Atractylodis Rhizoma (IcdP, DahA) were specifically tested at five times the LOQ. The recoveries and RSDs averaged 86.7 ± 15.1% and 8.5 ± 6.9%, respectively (Table S3). Therefore, this analytical method is considered to have good applicability based on good recoveries and precision shown in most samples, although the RSDs of BaP in Zingiberis Rhizoma and Aconiti Lateralis Radix Preparata were as high as 54.9% and 24.6%, respectively.

Monitoring of PAHs in Herbal Medicines
The newly developed GC-MS/MS method was applied to survey the pollution level of 8 PAHs in HMs purchased from the pharmaceutical markets in Korea. In this monitoring study, 459 samples of 51 different medicinal herbs were used. As a result, in only 8.3% (39 of 459 samples) and 37.3% (19 of 51 items), none of the 8 PAHs were detected above LOQs. Excluding samples in which all PAHs were not detected, the sum of 8 PAHs in individual 51 HMs ranged from 0.45 µg/kg in Anemarrhenae Rhizoma to 270.94 µg/kg in Zingiberis Rhizoma, with the average content ranging from 0.9 µg/kg in Araliae Continentalis Radix to 110.8 µg/kg in Coptidis Rhizoma Preparata cum Vinum. In terms of risk potential, it is worth mentioning that the sum of 8 PAHs exceeded the maximum allowable total amount of 4 PAHs of 50 µg/kg proposed by the EC [11] in 7.8% (36 samples) of the 459 samples and 37.3% (19 items) of the 51 kinds of HMs, although EFSA concluded that a system of eight substances (8 PAHs) would not provide much added value compared to a system of four substances (4 PAHs) [5]. In particular, 9 out of 19 items polluted by 8 PAHs above the LOQs belonged to the first group consisting of thermally processed herbal medicinal products (Table S4). Such PAHs might be inferred to be originated from the production process of those products. In addition, CHR was detected at 122.87 µg/kg in Aucklandiae Radix, which was the highest level among all the polluted levels of 8 PAHs.
In another study [39], a risk assessment was conducted utilizing information related to the average amount of 8 PAHs, dietary exposure and benchmark dose lower confidence limit (BMDL 10 ) based on the contamination level of the 8 PAHs in 51 HMs. HMs are usually available as decoctions or pills. Therefore, in that study, when consumed as a decoction, the average and highest conversion rates from ingredients to decoction were considered to be 6% and 20%, respectively. In the case of pills, a conversion rate of 100% was applied, assuming that the entire amount of pills was filled with HMs. As a result, by calculating the BMDL 10 (0.49 mg/kg body weight/day), daily exposure and dietary exposure (assuming the conversion rates of decoction and pill to be 20% and 100%, respectively), the margin of exposure (MOE) was found to be 10 5 or higher. According to the Food Standards Agency, United Kingdom [40], MOEs between 10 4 -10 6 were considered 'unlikely to be of concern', and those above 10 6 were considered 'highly unlikely to be of concern'. Therefore, although the potential harm of PAHs in HMs to Koreans could be considered low or negligible, action to minimize future exposure is still needed.

Chemicals and Reagents
Eight standards and two internal standards used in this study were purchased from Sigma-Aldrich (St. Louis, MO, USA). Another internal standard (BghiP-d12) was bought from Cambridge Isotope Laboratories (Wesel, Germany). Their abbreviations are shown in Table 1. HPLC-graded ethanol, hexane, dichloromethane and acetonitrile (ACN) were purchased from Merck KGaA (Darmstadt, Germany). N,N-dimethyl formamide of HPLC grade was manufactured by Daejung (Siheung, Korea). Sodium sulfate anhydrous was obtained from Junsei (Tokyo, Japan). LC-MS graded Water was purchased from Sigma-Aldrich (Oakville, ON, Canada). Hexane/dichloromethane mixture (3:1, v/v) was prepared by mixing 600 mL of hexane and 200 mL of dichloromethane. N,N-dimethyl formamidewater mixture (9:1, v/v) was prepared by mixing 100 mL of water with 900 mL of N,Ndimethyl formamide. A 1% sodium sulfate solution was prepared by weighing 10 g of sodium sulfate and dissolving it in 1 L of water.

Standard Solution Preparation
Three kinds of internal standards were used when calibrating 8 PAHs. CHR-d12 for BaA and CHR, BaP-d12 for BbF, BkF and BaP and BghiP-d12 for IcdP, DahA and BghiP were used as the internal standards. After weighing 2 mg of 8 types of standards and 3 types of internal standards individually, they were dissolved in 20 mL of ethanol to prepare 100 µg/mL of individual PAH standard stock solution and internal standard stock solution. The mixed standard solution of 8 PAHs was prepared at 1000 ng/mL and 500 ng/mL by taking a part of the volume of each standard stock solution and mixing them with the dichloromethane. The mixed internal standard solution was prepared at 1000 ng/mL by taking a part of the volume of each stock solution and mixing them with the dichloromethane. Through preliminary tests, the suitable calibration concentrations of 8 PAHs were designed into 1, 2, 5, 10, 25, 50 and 100 ng/mL. The calibration concentration of internal standard was also prepared to be 50 ng/mL. Calibration curves were established using 7 concentrations ranging from 1 to 100 ng/mL according to an internal calibration methodology.

Sample Preparation
Referring to the history of detection of BaP (unpublished data from MFDS), based on the high probability of generation of PAHs during thermal manufacturing process and the top 30% of production rates, 51 kinds of HMs were collected. Due to the efficiency of experiment and space constraints, HMs were divided into three groups. As the first group, 19   For the convenience of experimental process for validation, the samples were represented by 6 kinds of HMs according to appearance, matrix characteristics and stable chromatographic baseline seen in each group. Sanguisorbae Radix Carbonisatum of the first group, Achyranthis Radix, Scutellariae Radix and Angelicae Gigantis Radix of the second and Anemarrhenae Rhizoma and Asiasari Radix et Rhizoma of the last group were chosen as representatives ( Figure 5).
A total of 459 samples of 51 HMs, as described above, were purchased from herbal medicinal markets located in Daejeon, Korea, and used in this study. About 70% of samples were produced and imported from China. About 29% of samples were produced in Korea. Before carrying out any further experiments, the identities of the HMs were confirmed through the sensory test by the specialists. Voucher specimens were deposited at the National Institute of Food and Drug Safety Evaluation (Cheongju, Korea). The samples were chopped into small pieces and stored at −18 • C. About 500-600 g of each sample was equally divided into four parts to make them homogeneous. One of the sample parts was pulverized by a grinder (KSP-35, Koreamedi, Daegu, Korea).

372
A total of 459 samples of 51 HMs, as described above, were purchased from herbal 373 medicinal markets located in Daejeon, Korea, and used in this study. About 70% of sam-374 ples were produced and imported from China. About 29% of samples were produced in 375 Korea. Before carrying out any further experiments, the identities of the HMs were con-376 firmed through the sensory test by the specialists. Voucher specimens were deposited at 377 the National Institute of Food and Drug Safety Evaluation (Cheongju, Korea). The samples 378 were chopped into small pieces and stored at −18 °C. About 500-600 g of each sample was 379 equally divided into four parts to make them homogeneous. One of the sample parts was 380 pulverized by a grinder (KSP-35, Koreamedi, Daegu, Korea). 381

382
Liquid-liquid extraction was conducted in two ways depending on whether the sam-383 ples had essential oil or not. For the essential oil-free samples, N,N-dimethyl formamide-384 water mixture (9:1, v/v) was not added (hereinafter, referred to as non-DFA extraction 385 method). At the first step, 5.0 g of the ground sample was weighed in 250 mL bottle 386 (Nalgene, Thermo Fisher Scientific TM Waltham, MA, USA). After adding 100 mL of dis-387 tilled water thereto, sonication was performed for 60 min in an ultrasonic extractor (8510, 388 Branson, St. Louis, MO, USA) at room temperature. Then, 100 mL of hexane and 250 μL 389 of mixed internal standard solution (1000 ng/mL) were added thereto and were mixed 390 with a homogenizer (Mixer 2, OMNI, Georgia, USA) for 5 min. After sonicating again for 391 30 min, centrifugation was performed at 3220× g for 10 min in a centrifuge (Eppendorf,392 Hamburg, Germany). The supernatant, the hexane layer, was transferred to a separatory 393 funnel. As for the remains in the bottle, adding 50 mL of hexane, shaking for 5 min in a 394 shaker (Jisco, Seoul, Korea), centrifuging at 3220× g and transferring the supernatant to 395 the same separatory funnel were conducted twice in sequence. After combining all the 396 hexane layers obtained above in the separatory funnel, the hexane layer was washed by 397 adding 50 mL of water. After discarding the water, the hexane layer was dehydrated by 398 adding about 30 g of anhydrous sodium sulfate (Waco, Osaka, Japan) and filtered through 399 Whatman No. 41 filter paper (Whatman International Ltd., Buckinghamshire, UK). The 400

Sample Extraction
Liquid-liquid extraction was conducted in two ways depending on whether the samples had essential oil or not. For the essential oil-free samples, N,N-dimethyl formamidewater mixture (9:1, v/v) was not added (hereinafter, referred to as non-DFA extraction method). At the first step, 5.0 g of the ground sample was weighed in 250 mL bottle (Nalgene, Thermo Fisher Scientific TM , Waltham, MA, USA). After adding 100 mL of distilled water thereto, sonication was performed for 60 min in an ultrasonic extractor (8510, Branson, St. Louis, MO, USA) at room temperature. Then, 100 mL of hexane and 250 µL of mixed internal standard solution (1000 ng/mL) were added thereto and were mixed with a homogenizer (Mixer 2, OMNI, Atlanta, GA, USA) for 5 min. After sonicating again for 30 min, centrifugation was performed at 3220× g for 10 min in a centrifuge (Eppendorf, Hamburg, Germany). The supernatant, the hexane layer, was transferred to a separatory funnel. As for the remains in the bottle, adding 50 mL of hexane, shaking for 5 min in a shaker (Jisco, Seoul, Korea), centrifuging at 3220× g and transferring the supernatant to the same separatory funnel were conducted twice in sequence. After combining all the hexane layers obtained above in the separatory funnel, the hexane layer was washed by adding 50 mL of water. After discarding the water, the hexane layer was dehydrated by adding about 30 g of anhydrous sodium sulfate (Waco, Osaka, Japan) and filtered through Whatman No. 41 filter paper (Whatman International Ltd., Buckinghamshire, UK). The filtrate was concentrated to 2 mL under reduced pressure of 250 mbar at 45 • C in a water bath (Heidolph, Schwabach, Germany) ( Figure 2).
In the case of the samples containing essential oil, after transferring the supernatant hexane layer to a separatory funnel (F1), the step of adding N,N-dimethyl formamide-water mixture (9:1, v/v) to the separatory funnel was introduced (hereinafter, referred to as DFA extraction method). That is, 50 mL of N,N-dimethyl formamide-water mixture (9:1, v/v) was added to the hexane layer of F1, and after shaking it, the layer of the N,N-dimethyl formamide-water mixture was moved to another separatory funnel (F2). This step was repeated three times. After that, 100 mL of 1% sodium sulfate solution (Merck, Darmstadt, Germany) and 50 mL of hexane were put into F2. After shaking, the separated layer of hexane was transferred to the other separatory funnel (F3). Then, 35 mL of hexane was added to F2, shaken, and the hexane layer was also transferred to F3. This step, including adding 35 mL of hexane and transferring from F2 to F3, was repeated once more. Next, the hexane layer of F3 was washed by adding 50 mL of water, and the water was discarded. The remained hexane layer in the F3 was dehydrated by adding about 30 g of anhydrous sodium sulfate and filtered through Whatman No. 41 filter paper. The filtrate was concentrated to 2 mL under reduced pressure of 250 mbar at 45 • C in a water bath (Figure 2).

Sample Purification
In this study, SPE was applied for the purification procedure in basic. The 1g/6cc Florisil cartridge (Waters Oasis, Waters Corp., Milford, MA, USA) was activated by flowing out 10 mL of dichloromethane and 20 mL of hexane at a rate of 2-3 drops per second before use. After adding the concentrate to the activated cartridge, 20 mL of hexanedichloromethane mixture (3:1, v/v) was eluted at a rate of 2-3 drops per second. The eluate was blown away under nitrogen gas at 35 • C, and the residue was dissolved again with 5 mL of dichloromethane and filtered through 0.45 µm syringe filter for the instrumental analysis ( Figure 2).

GC-MS/MS Analysis
The GC-MS/MS system was coupled to Shimadzu Nexis GC-2030, Shimadzu GCMS-TQ8040NX tandem mass spectrometer in multiple reaction monitoring modes, Shimadzu AOC-20i Plus Auto-Injector and Shimadzu GCMSsolution workstation (Version 4.52) software (Shimadzu, Kyoto, Japan) in which Smart MRM (multiple reaction monitoring) creates MS/MS measurement program automatically. Analytes were separated with Agilent Technologies 30 m × 0.25 mm × 0.25 µm, DB-5MS column. Flow rate was 2 mL/min using helium as a carrier gas. The column oven temperature was programmed as follows: the initial temperature was 40 • C, increased to 220 • C at 30 • C/min, then ramped to 320 • C at 5 • C/min, which was held for 2 min. The total runtime was 28 min. The temperature of the injection port was maintained at 310 • C and injection volume was 1.0 µL in the splitless mode. The temperature of ion source and interface was 250 and 310 • C, respectively. Argon was used as a collision gas (2.0 mL/min).
The MRM for all scan transitions was performed in positive ion mode and at event times of 0.10 to 0.15 sec/cycle. The initial condition in terms of the m/z values, mainly referring to the study of Anderson et al. [34], was set by infusing reference solutions at a concentration of 1 mg/L in the scan mode of GC-MS/MS. Including the optimal precursor, ionic products, and collision energy for each PAH, MRM conditions were newly obtained by the infusion of reference solutions of the targeted compounds at a concentration of 100 ng/mL into the GC-MS/MS system. The MRM conditions for 8 PAHs established in this study appear in Table 1.

Statistical Evaluations
All experiments were carried out in triplicate unless stated otherwise. The results given are the average values. Statistical calculations and data analyses were performed using Microsoft Excel (version 2016, Microsoft, Redmond, WA, USA).

Matrix Effect
This study evaluated the matrix effect prior to conducting validation studies. When analysis using a mass spectrometer is performed, the ionization efficiency may be lowered or higher depending on the matrix of the blank sample [35]. Therefore, for each target sample, a standard solution containing a matrix by dissolving the standard in a blank sample and a standard solution in which the standard was dissolved in a solvent were prepared, respectively. Then, the degree of influence of the matrix was evaluated by comparing the slope of the calibration curve between the standard solution containing the matrix and the standard solution dissolved in the solvent. In addition, the blank sample solutions were prepared by the same pretreatment method described above. The calibration curves were achieved from seven-point concentrations ranging from 1 to 100 ng/mL. The matrix effect was evaluated by the following formula [37]: Matrix effect (ME, %) = (slope of the calibration curve in matrix/slope of the calibration curve in solvent) × 100.

Method Validation
Method validation followed the AOAC Official Methods of Analysis guidelines [38]. The validation study was designed by obtaining three different individual samples per matrix. Pulling out selectivity and specificity, linearity, limit of detection (LOD) and limit of quantification (LOQ), and accuracy and precision were performed to verify this study. For selectivity and specificity, the separation and discrimination of analytes in the presence of interfering substances were compared and confirmed by measuring the matrix blank, the PAH standard solution and the mixture of the standard and matrix. For linearity, a calibration curve was prepared with the concentration of the analyte and the peak area ratio by measuring the standard solution of 7 concentrations with an analysis device, and then the correlation coefficient was confirmed. LOD and LOQ were calculated based on the standard deviation of the response (σ) and the slope of the calibration curve (S) by the following equation: LOD = 3.3 × σ / S, LOQ = 10 × σ / S. For measuring the accuracy and precision, three different levels of standard solution were spiked into blank samples to achieve 2, 5 and 10 times (low, medium and high) the LOQ. After performing pretreatment with the spiked herbal medicine samples, the quantification analysis was conducted. The intra-day precision (repeatability, relative standard deviation [RSD]) and accuracy (recovery rates) were evaluated through the analysis on a single day. Along with these, the inter-day precision (reproducibility, RSDr) and accuracy were derived three times over three consecutive days (n = 9). To verify the validity of the analytical method, cross-validation between internal and external laboratories was conducted for selectivity, specificity, linearity, LOD, LOQ, accuracy and precision.

Monitoring of 8 PAHs in Herbal Medicine
Contamination levels of 8 PAHs in HMs were investigated by this modified GC-MS/MS analytical method. A total of 459 samples of 51 items (refer to Sample Preparation section) were collected, identified, stored, pretreated and analyzed for the polluted levels of 8 PAHs.

Conclusions
The determination method of 8 PAHs in six representatives of HMs was successfully improved by liquid and liquid extraction followed by florisil SPE purification and subsequent GC-MS/MS analysis. In particular, extraction steps were applied depending on whether the individual HMs contained essential oil components. With MRM conditions set for 8 PAHs, the MRM mode showed a high sensitivity for the 8 PAHs in the GC-MS/MS analysis. Meanwhile, by applying this method to monitor the contamination level of 8 PAHs in 459 samples of 51 kinds of HMs purchased from Korean pharmaceutical markets, any PAH above LOQs was detected in 91.7% of all samples. Based on the MOEs from these results, risks attributable to the PAHs in HMs could be interpreted as low or negligible, with action minimizing future exposure. Based on the good applicability, the developed analytical method is expected to be used as a standard analytical method not only for HMs but for agricultural products and foodstuffs that may contain PAHs.
Supplementary Materials: The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/molecules28093853/s1, Table S1: Matrix effects of 35 kinds of herbal medicines; Table S2: LOD, LOQ, recoveries, intra-day and inter-day RSDs and matrix effects (ME) obtained using the GC-MS/MS in Scutellariae Radix; Table S3: Recoveries and precisions of 45 kinds of herbal medicines for applicability excepting 6 representatives; Table S4: The item lists of the sum of 8 PAHs is over 50 µg/kg among 459 samples of 51 kinds of herbal medicines and maximum concentration in the same items.